Multiplicity correlations of intermediate-mass fragments with pions and fast protons 

in + Au 
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Low-energy vr"'' {E^ < 35 MeV) from '^^C+^^'^Au collisions at incident energies from 300 to 1800 
MeV per nucleon were detected with the Si-Si(Li)-CsI(Tl) calibration telescopes of the INDRA mul- 
tidetector. The inclusive angular distributions are approximately isotropic, consistent with multiple 
rescattering in the target spectator. The multiplicity correlations of the low-energy pions and of 
energetic protons {Ep > 150 MeV) with intermediate-mass fragments were determined from the 
measured coincidence data. The deduced correlation functions 1 -I- J? « 1.3 for inclusive event sam- 
ples reflect the strong correlations evident from the common impact-parameter dependence of the 
considered multiplicities. For narrow impact-parameter bins (based on charged-particle multiplic- 
ity), the correlation functions are close to unity and do not indicate strong additional correlations. 
Only for pions at high particle multiplicities (central collisions) a weak anticorrelation is observed, 
probably due to a limited competition between these emissions. Overall, the results are consistent 
with the equilibrium assumption made in statistical multifragmentation scenarios. Predictions ob- 
tained with intranuclear cascade models coupled to the Statistical Multifragmentation Model are in 
good agreement with the experimental data. 

PACS numbers: 25.70.Mn, 25.70.Pq, 25.75.Dw, 25.75.Gz 



I. INTRODUCTION 

Large transfers of energy are required in order to ini- 
tiate fragment emissions from excited heavy nuclei. The 
phase space for the socalled cracking mode, i.e. the dis- 
integration into several fragments, assumes its maximum 
for excitation energies of the order of 1 GeV for a ^ss^j 
nucleus 0, and the threshold for fragmentation is ex- 
pected near 3 MeV per nucleon 2] . The observation of a 
rise and subsequent fall of the fragment multiplicities Q 
in the decay of spectator residues produced at relativistic 
bombarding energies has demonstrated that energies of 
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this magnitude, and even exceeding it, can be reached 
iiili 

In these reactions at high energy, the excitation of A 
isobars represents an important mechanism for the trans- 
fer of energy from the relative motion of the colliding 
nuclei into other degrees of freedom. The strength of 
pion production with pion-to-participant ratios of up to 
10% at incident energies around 1 GeV per nucleon gives 
evidence for it 0. Pion production, at these energies, 
proceeds predominantly via two-step processes, with the 
excitation of a A resonance in a hard nucleon-nucleon 
collision being the first step. The excitation of A iso- 
bars is considered as important also for the heating of 
the spectator residues via reabsorption or multiple scat- 
tering of the A or pion As an illustration, 
the excitation energy of projectile spectators in ^^^Au -I- 



^^C collisions at 1 GeV per nucleon, as predicted by the 
Isabel Intranuclear-Cascade model, will be reduced by 
almost 20% if the A channels are artificially suppressed 
in the calculations 5]. Furthermore, data obtained with 
pion projectiles show that their rest energy is as effec- 
tive in transferring excitation e nerg y to the target as the 
same amount of kinetic energy jl4| . 

Direct experimental evidence for the mechanisms re- 
sponsible for the heating of the produced spectator sys- 
tems is difficult to obtain because of the high degree 
of equilibration that is apparent in their decay j^, Il5l | . 
Besides global parameters, such as the spectator mass 
and excitation energy, little measurable traces of the vi- 
olent initial reaction stages have been found in the decay 
patterns. Exceptions exist as, e.g., the sidewards en- 
hancement of the fragment emission in proton induced 
collisions with Ep > 10 GeV 0, 0|, which has been 
attributed to the kinematics of the primary nucleon- 
nucleon scatterings [l7l |. 

In this work, we analyze another signature which de- 
pends on the number and the character of the nucleon- 
nucleon interactions during the initial cascading process. 
We study the correlations of the multiplicities of charged 
pions and of energetic protons with the multiplicities of 
intermediate-mass fragments in reactions of ^^C -I- ^^^Au 
at incident energies from 300 to 1800 MeV per nucleon. 
The INDRA multidetector ,18| has been used, and the 
possibility to detect and identify low-energy 7r+ (iJ^r < 35 
MeV) and high-energy protons {Ep > 150 MeV) with the 
Si-Si(Li)-CsI(Tl) cahbration telescopes, in addition to 
the nearly An solid-angle coverage for light charged par- 
ticles and nuclear fragments, has been exploited. Pion- 
fragment and proton-fragment correlation functions were 
constructed from the measured singles and coincidence 
yields. 

The detected pions and protons indicate the occur- 
rence of hard primary collisions but, at the same time, 
carry away major amounts of energy. Pions and, in par- 
ticular, the low-energy pions considered here are believed 
to be emitted late, near the end of the primary reaction 
phase 0. The absorption effect of the spectators is evi- 
dent in the dependence of the relative pion rates on the 
mass of the colliding nuclei jl^. Therefore, in principle, 
both positive or negative correlations may be expected 
and are indeed observed or, at least, weakly indicated. 
As a main result of this work we find, however, that the 
deposited energy is equally spread over the degrees of 
freedom of the decaying spectator system, consistent with 
the equilibrium assumption made in statistical multifrag- 
mentation scenarios. Calculations performed with hybrid 
models , th e Liege nucleus-nucleus Intranuclear Cascade 
model T^l with percolation stage [20l l2ll | and the Dubna 
Intranuclear Cascade j23, , both coupled to the Sta- 
tistical Multifragmentation Model , are in good agree- 
ment with the experimental results. 
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II. EXPERIMENTAL DETAILS 

The experimental data were obtained as part of the 
INDRA experimental program carried out at the GSI 
Darmstadt which encompassed the investigations of the 
symmetric collision systems ^^^Au -I- ^^^Au ^A, 25, 2^ 
and Xe -I- Sn [g^ at intermediate energies and of the 
asymmetric systems "'^^C -|- -'^^^Au and ^^C -I- Sn at rela- 
tivistic bombarding energies. 

Beams of ^^C projectiles with incident energies 300, 
600, 1000 and 1800 MeV per nucleon provided by the 
heavy-ion synchrotron SIS were directed onto ^^^Au tar- 
gets with arcal densities of 2.0 or 2.8 mg/cm^, mounted in 
the center of the INDRA detection system. Annular veto 
detectors upstream and data sets measured with the tar- 
get being removed were employed in order to identify the 
interactions of the beam halo with the detectors and the 
mounting structures. The beam halo was negligible at 
300 MeV per nucleon but caused significant trigger rates 
at the highest incident energies, including events not ac- 
companied by a signal from one of the halo detectors 
which had outer diameters up to 10 cm. Corrections for 
this latter component were obtained from the no-target 
runs, normalized with respect to the number of registered 
halo events. 

The INDRA detector telescopes are arranged in 17 az- 
imuthally symmetric rings. Of these, the 8 more back- 
ward rings, covering the range of polar angles 45° < 
^'lab < 176° and consisting of ionization chambers fol- 
lowed by CsI(Tl) scintillators, are each equipped with a 
calibration telescope. These telescopes consist of pairs of 
a 80- /xm Si detector and a 2- mm Si (Li) detector which are 
mounted between the ionization chamber and the CsI(Tl) 
crystal of one of the modules of a ring 18]. They subtend 
a small solid angle Ail — 13.2 msr per telescope but pro- 
vide high resolution for the identification of charged par- 
ticles. The calibration of the 80-/im Si detectors has been 
obtained with a-particles from standard mixed-nuclide 
sources and from ^^^Pb sources providing the 6.1 MeV 
and 8.8 MeV hues from the ^^^Pb daughter decay. The 
calibration coefhcients of the subsequent detectors were 
then adjusted so as to obtain an optimum reproduction 
of the locations of identified particles in the measured 
AE — E maps. For this procedure, the energy loss and 
range tables of Ref. ^] were used. 

An example for the two-dimensional maps obtained 
with the signals recorded by the 2-mm Si(Li) detectors 
and the subsequent CsI(Tl) detectors is shown in Fig. ^ 
Besides the three strong groups of the hydrogen isotopes 
weaker groups corresponding to pions and to fast protons 
punching through the CsI(Tl) detector can be recognized. 
We assume that the branch of stopped pions contains pre- 
dominantly TT"*". The prompt from their decay add 
only a few MeV to the measured energy signal while the 
decay processes following the capture of negative pions 
by nuclei of the detector material may cause substantial 
additional energy deposits in the detector (cf. Ref. 29]). 
The branch of pions punching through the CsI(Tl) detec- 
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FIG. 1: Scatter plot of events in the plane of signals recorded 
with the 2-mm Si(Li) and the CsI(Tl) detectors of ring 15 
((fiab> = 134°) in the reaction ^^C + ^^''Au at E/A = 1000 
MeV. For clarity, the same spectrum is shown without (top) 
and with explanations (bottom). The three lines labelled p, 
d, t correspond to protons, deuterons and tritons stopped 
in the CsI(Tl) detector, respectively. The contours indicate 
the selected regions of stopped tt^ (full line) and of protons 
punching through the CsI(Tl) detector (dashed). 



tor contains charged pions of either polarity. It merges 
with the branch of punch-through protons and, therefore, 
can only be identified at fairly backward angles where the 
yields of fast protons are sufficiently weak. In the analy- 
sis, regions containing stopped 7r+ and fast protons were 
selected that excluded cross-contaminations. For ring 15 
{S\ah ~ 134°), this is illustrated in the bottom panel of 
Fig. n In the case of 7r+ , these regions extended from « 
6 MeV to between 30 and 40 MeV of energy deposited in 
the detectors of which 4 MeV were assigned to the pro- 
duced /i"*" . Slightly larger energy intervals were chosen at 
the more backward angles because of the reduced inter- 
ference from punch-through protons. The corresponding 
ranges of kinetic energies are about 2 MeV up to 35 MeV 
for 7r+ and about 150 MeV to 300 MeV for protons. For 
the subtraction of the background, caused mainly by re- 
actions in the detector material l30i| , equivalent regions in 
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the neighborhood of the particle branches were sampled. 

The observed total multiplicity of charged particles, 
Mc, has been chosen as a measure of the violence of the 
collision which, in first order, is a function of the im- 
pact parameter b. Following the geometrical prescrip- 
tion [31|, five impact-parameter bins were defined, each 
with a width Ab — 0.2 • 6o where 6o is the maximum im- 
pact parameter contributing at the trigger condition of 
at least three registered hits in the INDRA detection sys- 
tem. Note that neither the pions nor the punch-through 
protons are counted in Mc- 



III. EXPERIMENTAL RESULTS 

The mean multiplicities of intermediate-mass frag- 
ments with atomic number 3 < Z < 30, of stopped tt^, 
and of protons with 150 MeV < Ep < 300 MeV are 
shown in Fig. [3 as a function of the impact parameter 
derived from the total multiplicity of charged particles. 
All multiplicities show a general increase with decreasing 
impact parameter. The pion and proton multiplicities in- 
crease monotonically with the bombarding energy, most 
clearly visible in the central impact-parameter bins. The 
energy dependence is different for the fragment multi- 
plicities for which the highest values are reached at 600 
and 1000 MeV per nucleon. The maximum is expected 
for the most central collisions at these energies, accord- 
ing to the impact-parameter dependence observed in the 
reverse reaction '^^'^Au -I- '^^C [MlS^I- The measured val- 
ues of below 2, however, are only about half the maxi- 
mum multiplicities observed in the inverse-kinematics ex- 
periments which have no thresholds for the detection of 
projectile fragments H,I3- The main cause for this sub- 
stantial difference is the effective identification threshold 
for fragments of about 4 MeV per nucleon that was ap- 
plied in the present analysis, in addition to the effects 
of incomplete solid-angle coverage (« 90% of An) and 
of the finite thickness of the target and the shadowing it 
causes at angles near 9iah = 90° . Model calculations with 
parameters obtained from moving-source descriptions of 
the fragment emissions confirm that inefficiencies of this 
magnitude are to be expected for the present case of emis- 
sion from target spectators which are nearly at rest in the 
laboratory. These threshold effects depend sensitively on 
the small but finite momentum transfer to the decaying 
system and, thus, may vary with the incident energy. 
In addition, at 1800 MeV per nucleon and central colli- 
sions, the fragment multiplicity may already be past its 
maximum and in the regime of the fall of multifragment 
emission toward higher bombarding energies . 

The measured absolute pion multiplicities are rather 
small because of the small range of pion energies covered 
in the experiment. Extrapolated to a 47r solid angle, 
assuming isotropy (see below), the 7r+ multiplicity for 
-Ejr < 35 MeV reaches up to 0.35 for the most central 
bin at the highest bombarding energy. The multiplicity 
of energetic protons {Ep > 150 MeV) integrated over 




FIG. 2: Mean multiplicity of intermediate-mass fragments 
(top) and mean differential multiplicities dM/dSl of stopped 
tt"*" (middle) and of fast protons (bottom), as a function of the 
relative impact parameter b/bo for the reactions ^^C + ^^^Au 
at E/A ^ 300, 600, 1000, and 1800 MeV, distinguished by 
different symbols as indicated. The pion and proton multi- 
plicities are averaged over the angular range di^b > 45°. The 
data at E/A = 1800 MeV have been evaluated for impact 
parameters b < 0.6 bo. Only statistical errors are given which 
unless shown explicitly are smaller than the symbol size. 



FIG. 3: Angular distribution of the inclusive differential mul- 
tiplicity dM/dil of tt"*" for the interval of deposited energy E 
= 6-27 MeV for ^^C + ^^^Au at E/A = 600 MeV (triangles) 
and 1000 MeV (squares). The displayed errors are derived 
from the fluctuations of the background in the vicinity of the 
pion branches in the identification maps. The lines repre- 
sent the corresponding yields for E/A = 1000 MeV predicted 
by the Liege-Intranuclear-Cascade code (full line) and by a 
Maxwellian pion source with parameters T = 68 MeV and /3 
— 0.45 (dashed) which both fit well the vr^ spectra reported 



^'lab > 45° is of the order of 1. Their angular distribution 
is strongly forward peaked, however. 

The inclusive angular distributions of low-energy 7r+ 
for 600 and 1000 MeV per nucleon are shown in Fig. 
They were evaluated for a reduced interval of deposited 
energy, E = 6 - 27 MeV, which was found free of pro- 
ton contaminations in all spectra including those at the 
more forward angles. It corresponds to pion kinetic en- 
ergies £^7r ~ 2 - 23 MeV. Only the telescope belonging to 
ring 16 (6'iab « 149°) had a high AE threshold cutting 
into the pion branch and was not used for the angular 
distributions. The uncertainties displayed in this figure 
are those derived from the fluctuations of the background 
in the vicinity of the pion branches in the identification 
maps. These systematic errors are considerably higher at 
1000 MeV per nucleon because of the higher background 
level at this bombarding energy (note that only the sta- 
tistical errors are given for angle-integrated yields shown 
in other figures). 

The angular distributions at the two bombarding en- 
ergies are very similar, in first order isotropic, with a 
tendency for a slight rise toward the backward angles. 
This emission pattern is consistent with multiple rescat- 
tering and absorption processes in the target spectator 
which is practically at rest in the laboratory. It is in com- 
plete contrast to what an extrapolation of the data for 
high-energy pions would predict. Laue et ai, from their 
measurements for ^'^C + ^^"^Au at 1000 MeV per nucleon 
at 6'iab = 44° and 70° and with a threshold E,, > 150 
MeV, have identified a Maxwellian source with tempera- 
ture T = 68 MeV and velocity (3 = 0.45 which permits a 



good description of their spectra pjl|. Its extrapolation 
to backward angles and low energies leads to an angular 
distribution which is strongly forward peaked (Fig. |31 
dashed line). The effects of rescattering are clearly visi- 
ble in the angular distribution calculated with the Liege- 
Intranuclear-Cascade code (see Sect. IIV|) even though it 
is still slightly forward peaked (Fig. |3| full line). Overall, 
the predicted and measured differential tt"*" multiplicities 
are of comparable magnitudes. Extrapolated to An, the 
experimental inclusive 7r+ multiplicity for E^^ < 35 MeV 
is « 0.04 which amounts to about 10% of the inclusive 
and energy- integrated multiplicity reported by the KAOS 
collaboration |19l |. 

The multiplicity correlations between the detected low- 
energy tt"*" and intermediate-mass fragments, for 1000 
MeV per nucleon and after sorting into five impact- 
parameter bins, are shown in Fig. ^ The pion multi- 
plicity increases strongly with increasing centrality but 
exhibits no clearly recognizable additional dependence 
on the fragment multiplicity within a given impact- 
parameter bin. The correlation between these two quan- 
tities can, alternatively, be expressed in the form of the 
correlation function 



l + i? = 



(M^ ■ Mimf) 
(M^) ■ (Mimf) 



(1) 



constructed from the measured coincidences of pions and 
IMF's and from their single yields 33]. The brackets 
denote the mean values of the product and individual 
multiplicities obtained for a given data sample. For ori- 
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FIG. 4: Mean associated multiplicity of stopped Tr"*", averaged 
over Slab > 45°, as a function of the measured multiplicity of 
intermediate-mass fragments for the reaction ^'^C + ^^^Au 
at E/A = 1000 MeV. The data are sorted into 5 impact- 
parameter bins as indicated. 
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FIG. 5: Inclusive multiplicity correlation functions for 
stopped tt'*' and intermediate-mass fragments (top) and fast 
protons and intermediate-mass fragments (bottom) as a func- 
tion of the incident energy in ^'^C -I- ^^^Au. A rapidity y < 
0.3 was required for the intermediate-mass fragments. The 
displayed statistical errors are smaller than the data symbols 
with the exception of pions at E/A = 300 MeV. 



FIG. 6: Multiplicity correlation functions for stopped tt^ 
and intermediate-mass fragments as a function of the relative 
impact parameter b/bo for the reaction ^^C -|- ^^'^Au at E/A 
= 300, 600, 1000, and 1800 MeV. A rapidity y < 0.3 was 
required for the intermediate-mass fragments. The statistical 
errors are displayed. 
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FIG. 7: Multiplicity correlation functions for fast protons 
and intermediate-mass fragments as a function of the relative 
impact parameter b/bg for the reaction ^^C -I- ^^'^Au at E/A 
= 300, 600, 1000, and 1800 MeV. A rapidity y < 0.3 was 
required for the intermediate-mass fragments. The statistical 
errors are displayed. 



entation, the correlation function is unity (i? ~ 0) for 
uncorrelated emissions and, e.g., has a value 1 + R = 4/3 
for the special case of a strict proportionality of the two 
multiplicities and for a homogeneous population within 
finite multiplicity intervals that start from multipHcity 
zero. 

The correlation functions obtained for pions and frag- 
ments and for fast protons and fragments from the inclu- 
sive data sets at the four bombarding energies are shown 
in Fig. |S1 For the fragments, the conditions 3 < Z < 
30 and rapidity y < 0.3 were introduced so as to select 
intermediate-mass fragments emitted by the target spec- 



tator. The values of 1 + i? « 1.2 to 1.3 represent strong 
correlations of the considered multiphcities, as already 
evident from Figs. [21 E] The weak variations with the 
projectile energy are not systematic. No significant result 
is obtained for the pion-fragment correlation function at 
300 MeV per nucleon because the pion multiplicity is 
very small at this energy (Fig. [5J. 

The magnitude of the inclusive correlation functions, 
close to the value derived for the simple example given 
above, seems to indicate that they are dominated by the 
common variation of the considered multiplicities with 
the impact parameter. More pions and hard nucleon- 



nucleon scatterings are observed at more central impact 
parameters. If the contributions of these processes to 
the energy transfer to the spectator are completely equi- 
librated at the time of fragment emission, they will influ- 
ence both, the total charged particle multiplicity that is 
used for sorting and the multiplicity of fragments. This 
is what is observed. In order to test if additional event- 
wise correlations exist, the correlation functions were also 
evaluated for the finite impact-parameter bins generated 
from the total charged-particle multiplicity Mc (Figs. El 
O. For pions, the statistical uncertainties are large at 300 
MeV per nucleon and, at all energies, for the more pe- 
ripheral collisions for which the multiplicities are small. 
At 1800 MeV per nucleon, the beam halo introduces addi- 
tional uncertainties for the event samples at large impact 
parameter (6 > 0.6 6o). 

The exclusive correlation functions are, generally, 
much closer to unity, and most of the observed devia- 
tions from unity can be related to residual correlated de- 
pendences on impact parameter {M^) within the chosen 
finite intervals. At 1800 MeV per nucleon, e.g., the frag- 
ment multiplicity decreases with in the most central 
collisions but increases at0.4<6/feo<0.6 while the pion 
multiplicity increases monotonically with Mc. These neg- 
ative and positive correlations are reflected in the correla- 
tion function (Fig. El top left). Similarly, the overall ten- 
dency towards positive correlations at mid-peripheral im- 
pact parameters indicates that here, even within a given 
bin, the more violent collisions, characterized by a larger 
number of emitted pions and fast protons, imply also 
higher energy deposits and consequently more fragments 
for the produced spectator systems. These residual cor- 
relations should disappear, however, if the widths of the 
impact-parameter bins are reduced. As demonstrated in 
Fig. m this is also the case, seen most clearly, e.g., for 
protons in the interval 0.6 < 5/6o < 0.8. The results for 
the narrower bins are, on average, closer to unity than 
the value obtained by integrating over the full bin width. 

There are weak indications of anticorrelations between 
pions and fragments that might possibly result from a 
competition for excitation energy. At 1000 MeV per nu- 
cleon, e.g., the individual multiplicities increase with Ale 
(Fig. 121) which is not reflected in the values for the two 
most central bins. Choosing the narrowest possible bin- 
ning (open triangles in Fig. (HJ, a weighted average of 
1 + R = 0.96 ± 0.03 is obtained for b/ho < 0.4 from the 
eight bins covering this interval. This value is practically 
consistent with unity, i.e. no correlations, but may also 
be considered as giving a limit for possible anticorrela- 
tions caused by the competition of very slow pion and 
fragment emissions at the breakup stage. 

For energetic protons and fragments (Fig. O , the cen- 
trality dependence is the same at all four energies and, 
overall, very similar to that observed for the pions (Fig. 
|S1). The correlation functions are unity for central colli- 
sions, 1 + R = 1.00 ± 0.01 for 6 < 0.2 bo and averaged 
over all incident energies. With increasing impact pa- 
rameter, they rise monotonically up to maximum values 
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FIG. 8: Multiplicity correlation functions for stopped tt^ 
(top) and fast protons (bottom) with intermediate-mass frag- 
ments for the reaction ^^C -I- ^^^Au at E/A = 1000 MeV and 
for impact-parameter bins of widths Ab/bg — 0.05 (open tri- 
angles), 0.1 (open squares), and 0.2 (dots, same as in Figs. |S1 
ITJ. A rapidity y < 0.3 was required for the intermediate-mass 
fragments. The statistical errors are displayed. 



of 1-1- i? w 1.4 but then drop in the bin of largest impact 
parameter to a mean value 1+R w 1, indicating uncorre- 
lated emissions of the two species in the most peripheral 
reactions. This loss of correlation seems to be present 
also in the pion case but is seen less clearly because of 
the very low emission rates. 



IV. DISCUSSION 

Two hybrid models were employed for the comparison 
of the measured multiplicity correlations with theoreti- 
cal expectations. The Liege relativistic nucleus-nucleus 
cascade |0| has recently been coupled to a percolation 
stage used for identifiying excited clusters and nuclear 
fragments in the nucleon distribution after the cascading 
stages of the reaction j23i 1^ ■ Their subsequent decay 
is then followed with the Statistical Multifragmentation 
Model (SMM) in the version described in Ref. 0. In 
the Liege-Intranuclear-Cascade code, the pions are pro- 
duced through the reactions NN ^ NA and A ^ TriV. 
The needed cross-sections have been fitted to experimen- 
tal data when available Only the A(1232) reso- 
nance is considered. The cross-section aNA^NN is mul- 
tiplied by a factor of 3 from the one deduced from the 
detailed balance of the reaction NN — > NA as discussed 
in 13 [s^. A test of the code has been performed by 
calculating the energy spectra of tt"*" emission for the ^^C 
+ ^^^Au reaction at 1000 MeV per nucleon and for the 
angles 6'iab = 44° and 70° chosen in the KAOS experi- 
ment [Tgj . Nuclear surfaces are sharp in the version used 
here, and the absolute reaction cross-section has to be 
estimated. With a normalization with respect to a to- 
tal reaction cross-section of 3 barns [l^, the obtained 
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FIG. 9: Predictions for tt^ kinetic-energy spectra at Siab ~ 
44° and 70° for ^^C + ^^''Au at E/A = 1000 MeV obtained 
with the Liege-Intranuclear-Cascade (histograms) in compar- 
ison with the data from the KAOS experiment (squares, Ref. 
figll'). The calculations are normalized assuming a total reac- 
tion cross section of 3 b. 

agreement is very satisfactory (Fig.EI). 

In a second approach, the intranuclear cascade (INC) 
model developed in Dubna was used for the simulation 
of the initial stage of the collision "22*, . The INC de- 
scribes the process of the hadron-nuclcon collisions inside 
the target nucleus. High energy products of these inter- 
actions, including pions, are allowed to escape while low 
energy products are assumed to be trapped by the nu- 
clear potential of the target system. At the end of the 
cascade, a residue with a certain mass, charge and exci- 
tation energy remains which then can be used as input 
for the statistical description of the fragment production. 
For this second stage also the Statistical Multifragmen- 
tation Model has been used. The model provides the 
option of inserting a preequilibrium stage between the 
INC and fragmentation stages 36.j . 

The predictions of both models are in good agreement 
with the data (Fig. HUJl. There is no significant difference 
caused by the fact that the intermediate configuration of 
the system is characterized in different ways in the two 
approaches, by the actual phase space configurations of 
nucleons in the Liege/percolation model while only the 
mass and excitation energy of a single spectator source 
are transferred to the SMM in the Dubna/SMM case. 
Only for the most peripheral collisions, large correlations 
are predicted by the models but not observed in the ex- 
periment. 

The correlations of the mean excitation energy with the 
total number of pions or fast protons, as obtained with 
the Dubna INC model, arc shown in Fig. [ill The INC 
output was sorted into bins of 2- MeV width in excitation 
energy, and the mean values were determined individu- 
ally for the four bins spanning the range from to 8 MeV 
per nucleon. If an intermediate stage of pre-equilibrium 
emission subsequent to the cascading stage is added the 
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FIG. 10: Predictions for ^^C ^ ^^"^ Ku at EjA = 1000 MeV 
obtained with the Liege-Intranuclear-Cascade with percola- 
tion stage (full line) and the Dubna-Intranuclear-Cascade 
(dashed line), both coupled to the Statistical Multifragmenta- 
tion Model, in comparison with the experimental multiplicity 
correlations of pions (top) and fast protons (bottom) with 
intermediate-mass fragments (squares). The predictions de- 
rived from the smooth interpolation described in the text are 
given by the dotted line. 
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FIG. 11: Results of calculations with the Dubna INC model, 
without (left panel) and with preequilibrium stage (right 
panel), for ^^C -|- "'^Au collisions at EjA = 1000 MeV: mean 
value of the excitation energy of the spectator residue, evalu- 
ated within bins of width 2 MeV per nucleon, as a function of 
the multiplicities of pions (full line) and fast protons (dashed). 



remaining excitation energy will be reduced. The corre- 
lations obtained in this case are shown in the right panel 
of Fig. ^2 These mean excitation energies are above or 
below the middle of the bin if the distribution of events 
is not homogeneous and skewed either on the low or on 
the high energy side. Their correlation with the mul- 
tiplicity of pions or fast protons thus indicates whether 
the observation of these particles can have an additional 
influence on the distribution of excitation energies once 



a bin has been selected, e.g., on the basis of the total 
charged-particle multiplicity. The obtained correlations 
are positive at low excitation energy and gradually dis- 
appear as the excitation energy is raised. In the range of 
low excitation energies close to the threshold for multi- 
fragment emissions, there is also a positive correlation 
between the multiplicity of intermediate-mass fragments 
and the excitation energy, so that the expected multiplic- 
ity correlations are positive, as predicted by the models 
(Fig. I1U|) . At higher excitation energies, the correlations 
are flat indicating that the multiplicity of the considered 
particles does not represent an additional weighting fac- 
tor for the distribution of excitation energies. 

The calculations suggest that the variations of the de- 
posited excitation energy within the chosen finite inter- 
vals of charged-particle multiplicity (impact parameter) 
are sufficient to explain the gross behavior of the mea- 
sured correlation functions. This conclusion has been 
further tested with model calculations solely based on 
a phenomenological analysis of the experimental multi- 
plicities. For this purpose, smoothed interpolations of 
the measured fragment, pion and fast-proton multiplici- 
ties, as functions of the charged-particle multiplicity Mc, 
were generated. The product multiplicities were approx- 
imated by the products of the corresponding singles mul- 
tiplicities. These quantities were then integrated over 
the ranges of Mc that correspond to a given impact- 
parameter bin, with weights given by the cross sections 
associated with the individual values of Mc within that 
bin. Correlation functions were finally calculated from 
the integrated singles and coincidence multiplicities (dot- 
ted fines in Fig. [Tn|l . 

In this approximation, the correlations will trivially 
disappear {R = 0) if the bin width is reduced to a single 
unit of Mc- In bins of finite width, deviations from unity 
may result from correlated Mc dependences of the sin- 
gles multiplicities of the considered particle species. Ad- 
ditional event-wise correlations in the experimental data 
would then appear as deviations from the so-obtained 
predictions. This, however, does not seem to be the case. 
The overall trend of the data is rather well reproduced, 
with the exception of the most peripheral bin (Fig. 1101) . 

The anticorrelation for pions in the most central event 
sample may indicate a weak competition between slow 
pions and fragments, on a level not larger than 10%. It 
may also exist at larger impact parameters where it is 
not visible in the presence of the residual positive corre- 
lation caused by the finite bin widths. The loss of cor- 
relations in very peripheral collisions is not predicted by 
the calculations. This narrows the class of possible expla- 
nations down to processes not included in the INC model 
as, e.g., scattering from correlated structures in the tar- 
get nucleus |37j . The cascade processes produce fast 
protons with high associated transfers of momentum as 
a result of multiple scatterings from individual nucleons 
inside a nucleus. The excitation energy of the produced 
residue is proportional to the number of such interac- 
tions and, therefore, is expected to be high in the case 
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of energetic protons at backward angles. The scattering 
from correlated nucleons (clusters) can produce large mo- 
mentum transfers with relatively low transfer of energy, 
corresponding to the high effective mass of such processes 
i37i | . This direct process, even if relatively rare, may con- 
tribute significantly to processes with high momentum 
transfers in very peripheral collisions (small multiplici- 
ties). Taking into account, however, that the multiplic- 
ities are extremely small and the systematic effects due 
to halo events largest in this impact-parameter bin, any 
such interpretation has to be considered speculative at 
this time. 



V. CONCLUSION 

The investigated correlations connect the production 
of particles considered as representative for the primary 
cascading stage of the reaction with the intermediate- 
mass fragments emitted in the spectator decay in ^^C -I- 
i^'^Au reactions at 300, 600, 1000, and 1800 MeV per 
nucleon. The strong correlations obtained for inclusive 
data sets reflect the common impact-parameter depen- 
dence of the considered multiplicities. More violent colli- 
sions, characterized by a larger number of emitted pions 
and fast protons, also imply higher energy deposits for 
the produced spectator systems. 

No strong evidence for additional direct correlations 
has appeared from the exclusive correlations functions, 
studied in their dependence on the multiplicity of charged 
particles. A weak anticorrelation between pions and frag- 
ments, limited to a few percent, may be indicated in cen- 
tral collisions. It is most likely due to a weak competi- 
tion between these emissions, caused by the about 160 
MeV of energy that is lost from the system with a slow 
pion. This is not observed for the fast protons which 
carry away equal amounts of energy but, apparently, are 
emitted earlier. 

The situation is less clear for very peripheral collisions 
for which a tendency toward uncorrelated emissions of 
intermediate-mass fragments and either slow pions or fast 
protons is observed, contrary to the model predictions. 
The scattering of nucleons from cluster structures has 
been given as an example for a mechanism associated 
with high momentum transfers and small energy transfers 
that is not considered in the models. 

Overall, the observed multiplicity correlations confirm 
the equilibrated nature of the produced spectator con- 
figurations prior to their decay. The multiplicity corre- 
lations of intermediate-mass fragments with either pions 
or fast protons are implicitly given by their correlations 
with the charged-particle multiplicity. This conclusion is 
supported by calculations based on a smooth interpola- 
tion of the fragment, pion and fast proton multiplicities 
as a function of the charged-particle multiplicity. 

Two-stage hybrid models reproduce the observations 
quantitatively, independently of their differences in the 
mode of coupling the two reaction stages. In the 



INC/SMM model, the parameters transferred to the sec- 
ond stage are restricted to the global variables of mass 
and excitation energy of the produced spectator system, 
and direct correlations are excluded since pions or very 
energetic protons are not produced in the second stage. 
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